Exact edge singularities and dynamical correlations in spin- 1/2 chains 
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Exact formulas for the singularities of the dynamical structure factor, S zz (q, oj), of the S = 1/2 xxz 
spin chain at all q and any anisotropy and magnetic held in the critical regime are derived, expressing 
the exponents in terms of the phase shifts which are known exactly from the Bethe ansatz solution. 
We also study the long time asymptotics of the self-correlation function (0\Sj(t)Sj(G)\0). Utilizing 
these results to supplement very accurate time-dependent Density Matrix Renormalization Group 
(DMRG) for short to moderate times, we calculate S zz (q,uj) to very high precision. 
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The "xxz" S =1/2 spin chain, with Hamiltonian 
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is one of the most studied models of strongly correlated 
systems. It is equivalent by a Jordan- Wigner transforma- 
tion to a model of interacting spinless fermions, with the 
corresponding Fermi momentum kp — vr(l/2 + (0|S'J|0)) 
1] . The model with A = 1 describes Heisenberg antifer- 
romagnets. The regime < A < 1 is also of experimental 
interest; for example, the model with A = 1/2 can be re- 
alized in 5=1/2 spin ladders near the critical field Q. 
In optical lattices, it should be even possible to tune the 
anisotropy A and explore the entire critical regime [3j]. 

While some aspects of the model have been solved for 
exactly by Bethe ansatz it has been very difficult to 
obtain correlation functions that way. Field theory (FT) 
methods give the low energy behavior at wave-vectors 
near and 2kp [H ■ From the experimental viewpoint [j| , 
a relevant quantity is the dynamical structure factor 
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This is the Fourier transform of the density correlation 
function in the fermionic model. For A = 1 and h = 0, 
the exact two-spinon contribution to S zz (q,u>) was ob- 
tained from the Bethe ansatz [6( , partially agreeing with 
the Miiller conjecture Q. More recently a number of 
new methods have emerged which now make this prob- 
lem much more accessible. These include time-dependent 
DMRG [8L Isl, lOj , calculation of form factors from Bethe 
ansatz [ll|,|l2| and new field theory approaches which go 
beyond the Luttinger model lj, [l4j . The results point 
to a very nontrivial line shape at zero temperature for 
S zz (q,u>) of the xxz model 14 1 and of one-dimensional 
models in general [l3| . In the weak coupling limit A< 1 
and for small q, the singularities at the thresholds of the 
two-particle continuum have been explained by analogy 
with the x-ray edge singularity in metals (l3| . 



In this Letter we combine the methods of Ref. [13[ 
with the Bethe ansatz to investigate the singularity expo- 
nents of S zz (q,uj) for the xxz model for finite interaction 
strength A and general momentum q. In addition, we de- 
termine the exponents of the long-time asymptotics of the 
spin self-correlation function, which is not dominated by 
low energy excitations. We check our predictions against 
high accuracy numerical results calculated by DMRG. 

In the non-interacting, A = case, only excited states 
with a single particle- hole pair contribute to S zz {q,ui). 
All the spectral weight is confined between the lower and 
upper thresholds u>l,u(q) of the two-particle continuum. 
The choices of momenta corresponding to the thresholds 
depend on both kp and q. For zero field, kp = tt/2, 
uoL.{q) for any q > is defined by the excitation with a 
hole at k\ = ir/2 — q and a particle right at the Fermi 
surface (or a hole at the Fermi surface and a particle at 
&2 = 7r/2 + q), while ujjj is defined by the symmetric 
excitation with a hole at k\ = tt/2 — q/2 and a particle 
at ki = 7r/2 + q/2. For finite field and q < \2kp — 7r|, 
^L,u{q) arc defined by excitations with either a hole at 
kp and a particle at kp + q or a hole at kp — q and a 
particle at kp. For h ^ and q > \2kp — tt|, there is even 
a third "threshold" between wr and tujj where S zz (q,ui) 
has a step discontinuity (see [3]). 

For A / 0, S zz (q,ui) exhibits a tail associated with 
multiple particle-hole excitations 14j. However, the 



thresholds of the two-particle continuum are expected 
to remain as special points at which power-law singular- 
ities develop [13|. In order to describe the interaction 
of the high energy particle and/or hole with the Fermi 
surface modes, we integrate out all Fourier modes of the 
fermion field tp(x) except those near ±kp and near the 
momentum of the hole, k\, or particle, k%, writing 
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Linearizing the dispersion relation about ±kp we obtain 
rclativistic fermion fields which we bosonize in the usual 
way [l|. We also expand the dispersion of the d 1)2 parti- 
cles around k = k%% up to quadratic terms. This yields 
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the effective Hamiltonian density 
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This Hamiltonian describes a Luttinger liquid coupled to 
one or two mobile impurities (TBI . [Tq . In the derivation 
of Eq. (|4]) from Eq. (fTJ), we drop terms of the form 
(dad a ) 2 because we only consider processes involving a 
single d\ and/or a single d,2 particle. Here <pR t L are the 
right and left components of the rescaled bosonic field. 
The long wavelength fluctuation part of S z is given by 

Sj ~ y/K/2n (d x (fR + d x ipi). The spin velocity v and 
Luttinger parameter K are known exactly from the Bethe 
ansatz 4J- For zero field, v = (it /2)\/l — A 2 / arccos A 
and K = [2 - 2 arccos(A)/7r)] _1 (we set J = 1). To 
first order in A, the coupling constants describing the 
scattering between the d particles and the bosons are 
k r l — 2A[1 — cos(kp^fk a )]- The direct di-d-z interaction 
V\2 is also of order A. The exact values of Krl play a 
crucial role in the singularities and will be determined 
below. 

We may eliminate the interaction between the d parti- 
cles and the bosonic modes by a unitary transformation 
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with parameters 7^ L = k r l /(v =p u a ). In the resulting 
Hamiltonian H = WHU, <pr.l are free up to irrelevant 
interaction terms [TB| . As in the x-ray edge problem, 7^ L 
may be related to the phase shifts at the Fermi points due 
to the creation of the high energy d a particle. 

Fortunately, we have access to the high energy spec- 
trum of the xxz model by means of the Bethe Ansatz. 
Following the formalism of Ref . 4j| , we calculate the fi- 
nite size spectrum from the Bethe ansatz equations with 
an impurity term corresponding to removing (adding) 
a particle with dressed momentum k\ = k(X%) (k 2 = 
fc(A2)), where Ai^ are the corresponding rapidities. The 
term of 0(1) yields e a = e(k a ), the dressed energy of 
the particle. For zero field, we have the explicit formula 
e(k) = —vcosk. The excitation spectrum for a single 
impurity to 0(1/ L) reads 
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with a conventional notation for AN, D and n± 
phase shifts nf mp and df mp are given by 



d^ pZnpW> 



(6) 
The 

(7) 



imp 



d\ 



dX tmpK . 



(8) 



where B is the Fermi boundary and pf mp (X) is the solu- 
tion to the integral equation 
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where 6(A) = ilog[sinh(iC + A)/sinh(i£ — A)l, with 
A = — cos£, is the two-particle scattering phase [J], and 
$ 1,2 (A) = TdO(X - X 1<2 )/dX. The spectrum of Eq. © 
describes a shifted c = 1 conformal field theory (CFT). 
The scaling dimensions of the various operators can then 
be expressed in terms of K, nf mp and df mp . In the ef- 
fective model (j4]), the shift is introduced by the unitary 
transformation of Eq. ([5]), which changes the boundary 
conditions of the bosonic fields. The equivalence of the 
two approaches allows us to identify 



iRLh 



77 a +2Kd a 

v imp ^imp m 



(10) 



The phase shifts can be determined analytically for zero 
magnetic field. In this case, B — 
0. Moreover, by integrating Eq. 



00 and we have df mp 
(HU) over A we find 
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Once the exact phase shifts are known, the exponent 
for the (lower or upper) threshold determined by a sin- 
gle high energy particle can be calculated straightfor- 
wardly. For example, for a lower threshold defined by 
a deep hole, = — e(^F — <?), the correlation func- 

tion (d\ip R (t, a^Jjd-L (0, 0)) can be factorized into a free 
d\ propagator and correlations of exponentials of ifRL- 
After Fourier transforming, we find that near the lower 
edge S zz (q,ui) ~ [w — uji J (q)]~ >1 with exponent [17 1 
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For h — > 0, we use Eq. pT|) and obtain 
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independent of the momentum of the hole. This form for 
the lower edge exponent had been conjectured long ago 
by Miiller et al. [7[. It agrees (up to logarithmic correc- 
tions) with the exponent of the two-spinon contribution 
to S zz (q,ui) for the Heisenberg point (K = 1/2) [6|. 

The general result of Eq. (fl2l is consistent with the 
weak coupling expression for /i [13j. To first order in A, 
Eq. (fT2|) reduces to 
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For kp ^ 7r/2, we expand for q -C kp and get fi ~ 
mAq/n, where m = (cos kp)" 1 (c.f. (lj|). For kp = tt/2, 
we obtain fi w 2A/ir, which is 1 — K to O(A). Note the 



3 



cancellation of the q dependence of and v — u\ in the 
latter case. Momentum-independent exponents have also 
been derived for the Calogero-Sutherland model [HI . 

We now consider a threshold defined by high-energy 
particle and hole at fci.2 = tt/2 =p q/2. The relevant cor- 
relation function is the propagator of the transformed 
d\d x . For simplicity, here we focus on the zero field 
case, in which £2 = — £1 = u sin(g/2), ui = u\ and 
[v sin(g/2)] _1 . Particle-hole symmetry 
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then implies that 7^ L = 7^ L and d\d 1 is invariant 
under the unitary transformation of Eq. ([5]). In the 
nonintcracting case, there is a square root singularity 
at the upper threshold due to the divergence of the 
joint density of states: S zz (q,to) oc \J mi / [ujjj {q) — w\ for 
u) ps u>u(q) = 2w sin(q/2) [7|. For A 7^ 0, we need to treat 
the direct interaction V12 between the particle and the 
hole, which is not modified by U . This problem is anal- 
ogous to the effect of Wannier excitons on the optical 
absorption rate of semiconductors [3, |2(J. This sim- 
ple two-body problem can be solved exactly for a delta 
function interaction. The result is that the upper edge 
exponent changes discontinuously for A 7^ 0: the square 
root divergence turns into a universal (for any q and A) 
square root cusp, S zz (q,u>) oc i/wyfa) — u>. This behav- 
ior contradicts the Muller ansatz [7|, but is consistent 
with the analytic two-spinon result for A = 1 6] . Unlike 
the original exciton problem, a bound state only appears 
for Via < (A < 0) [2l|, because the particle and hole 
have a negative effective mass. For the upper edge 

cusp should intersect a high-frequency tail dominated by 
four-spinon excitations as proposed in [22| ■ This picture 
must be modified for h 7^ 0, since then 7^, L 7^ Jrl and 
one needs to include the bosonic exponentials. The upper 
edge singularity then becomes A- and ^-dependent. The 
eneral finite field case, including the middle singularity 
for q > \2kF — tt|, will be discussed elsewhere. 
We can apply the Hamiltonian of Eq. ^ to study the 
self-correlation function G(t) = (0|5|(t)5|(0)|0). Even 
in the nonintcracting case, the long time asymptotics is 
a high energy property, since it is dominated by a saddle 
point contribution with a hole at the bottom and a par- 
ticle at the top of the band (23)]. In this case, ki — and 
k 2 = 7r and d\^ p vanish by symmetry (7^ = 7£). Here 
we restrict to zero field, but the method can be easily 
generalized. For h — and A > 0, G(t) takes the form 
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where W = — e(0) = v. The last two terms are the 
standard low-energy contributions, with a = 2K. The 
amplitudes B3 and B4 are known 24|. The first term is 
the contribution from the hole at the bottom of the band 
and the particle at fcp = 7r/2, with exponent 

t? = (1 + K)/2 + (1 - n\ mp f/2K = K + 1/2. (16) 



The term oscillating at 2W comes from a hole at k = 
and a particle at k = ir. For A = 0, we have 
772 = I- The exponent 772 is connected with the sin- 
gularity at the upper threshold of S zz (q,ui) by G(t) ~ 
J dioe lwt J dqS zz (q,uj) for q ps 7r and w w u>u{^) = 2u. 
Due to the discontinuity of the exponent at wp, 772 jumps 
from 772 = 1 to 772 = 2 for any nonzero A. This behavior 
should be observed for t ^> l/(miVi 2 ) ~ 1/A 2 . As a re- 
sult, the asymptotics of G{t) is governed by the exponent 
j] < 3/2 for < A < 1. For A < 0, we must add to Eq. 
(TTS)) the contribution from the bound state. 

We can also study S zz {q 1 uo) with time-dependent 
DMRG (tDMRG) The tDMRG methods di- 

rectly produce S zz {x, t) and its spatial Fourier transform 
S zz (q,t) for short to moderate times. This information 
nicely complements the asymptotic information available 
analytically. The DMRG calculation begins with the 
standard finite system calculation of the ground state 
4>(t = 0) on a finite lattice of typical length L — 200-400, 
where a few hundred states are kept for a truncation er- 
ror less than 10 -10 . One of the sites at the center of the 
lattice is selected as the origin, and the operator Sq is 
applied to the ground state to obtain a state t/j(t = 0). 
Subsequently, the time evolution operator for a time step 
r, exp(i(H — Eq)t) where Eq is the ground state energy, 
is applied via a fourth order Trotter decomposition [nj to 
evolve both <j>{t) and ip(t). At each DMRG step centered 
on site j we obtain a data point for the Green's function 
G(t,j) by evaluating {(^{t^S^iplt)) . As the time evolu- 
tion progresses, the truncation error accumulates. The 
integrated truncation error provides a useful estimate of 
the error, and so longer times require smaller trunca- 
tion errors at each step, attained by increasing the num- 
ber of states kept m. The truncation error grows with 
time for fixed m, and is largest near the center where 
the spin operator was applied. We specify the desired 
truncation error at each step and choose m to achieve 
it, within a specified range. Typically for later times we 
have m ps 1000. Finite size effects are small for times 
less than (L/2)/v. We are able to obtain very accurate 
results for G(t,j), with errors between 10 -4 and 10~ 5 , 
for times up to Jt ~ 30-60. 

For Jt > 10 - 20, we find the behavior of S zz (q, t) and 
G(t) is well approximated by asymptotic expressions, de- 
termined by the singular features of S zz (q, to) and G(w). 
By utilizing the leading and subleading terms for each 
singularity, we have been able to fit with a typical error in 
S zz (q,t) or G(t) for Jt - 20-30 between 10" 4 and 10~ 5 . 
We can fit with the decay exponents determined analyt- 
ically or as free parameters to check the analytic expres- 
sions. Table U shows the comparison between the expo- 
nents for G(t) extracted independently from the DMRG 
data and the FT predictions. In all cases the agreement is 
very good. By smoothly transitioning from the tDMRG 
data to the fit as t increases, we obtain accurate results 
for all times. A straightforward time Fourier transform 
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TABLE I: Exponents for the spin self-correlation function 
G(t) for h — 0. The parameters W, n, rj2 and a were ob- 
tained numerically by fitting the DMRG data according to 
Eq. (I15|l . These are compared with the corresponding FT 
predictions (with v and K taken from the Bethe ansatz). 
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the particle-hole symmetric zero field case, we showed 
that the lower edge exponent is g-independent and the 
( "exciton-like" ) upper edge has a universal square root 
singularity. The combination of analytic methods with 
the tDMRG overcomes the finite t limitation on the res- 
olution of the tDMRG and can be used to study dynamics 
of other one-dimensional systems (integrable or not). 
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Glazman, N. Kawakami, M. Pustilnik and J. Sirker for 
discussions. We acknowledge the support of the CNPq 
grant 200612/2004-2 (RGP), NSERC (RGP, IA), NSF 
under grant DMR-0605444 (SRW) and ClfAR (IA). 



with a very long time window yields very accurate high 
resolution spectra. Examples of line shapes obtained this 
way are shown in Fig. [TJ We also did DMRG for the hole 
Green's function for the fermionic model corresponding 
to Eq. JT]) , obtaining good agreement with the predicted 
singularities from the x-ray edge picture. 

We have not seen any exponential damping of the 772 
term in G(t) for A > 0. This suggests that the singularity 
at the upper edge is not smoothed out in the integrable 
xxz model, even when the stability of the excitation is 
not guaranteed by kinematic constraints [2|| . Integrabil- 
ity also protects the singularity at u>u for finite field, as 
implied by the CFT form of the spectrum in Eq. ((6|) . 

In conclusion, we presented a method to calculate the 
singularities of S zz (q,u>) for the xxz model. The ex- 
ponents for general anisotropy, magnetic field and mo- 
mentum can be obtained by solving the Bethe ansatz 
equations which determine the exact phase shifts. For 
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FIG. 1: (Color online). DMRG results for S zz (q,u) versus 
uj for q = 7r/2, h — and several values of anisotropy A. 
The line shapes for A > show a divergent x-ray type lower 
edge and a universal square-root cusp at the upper edge. The 
curve for A < shows a bound state above the upper edge. 
The width of the peak is very small for small |A|. 
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